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Abstract 

Bolgiano scaling in Boussinesq turbulence is analyzed using the Yakhot-Orszag renormal- 
ization group. For this purpose, an isotropic model is introduced. Scaling exponents axe 
calculated by forcing the temperature equation so that the temperature variance flux is 
constant in the inertial range. Universal amplitudes associated with the scaling laws are 
computed by expanding about a logarithmic theory. Connections between this formalism 
and the direct interaction approximation are discussed. It is suggested that the Yakhot- 
Orszag theory yields a lowest order approximate solution of a regularized direct interaction 
approximation which can be corrected by a simple iterative procedure. 


I. Introduction 

Application of the Kolmogorov scaling theory of hydrodynamic turbulence to Boussi- 

nesq turbulence suggests two regimes. The simpler possibility is that the gravitational 
coupling can be neglected. Then the temperature is a passive scalar and the velocity and 
temperature fluctuations have Kolmogorov spectra 

E ~ e I. 2 / 3 * jb -5 / 3 
E t ~ Ne~ 1/3 k~ 5 ^ 

where e is the kinetic energy dissipation rate, Et is the spectrum for mean square temper- 
ature fluctuations, and N is the corresponding dissipation rate. In the passive regime, the 
correlation < u\T' > vanishes by symmetry, and the corresponding spectrum Eh is zero. 
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Bolgiano 1 identified the second possibility, that velocity and temperature fluctuations 
are deter min ed by N and g. Then dimensional analysis leads to 

E ~ g 4/5 jV 2/5 Ar 11/5 

EH-g^N^Sk-*/ 5 

E T ~g- 2 ' 5 N 4 ' 5 k - 7 / 5 

The velocity spectrum falls off more quickly at small scales than a Kolmogorov spectrum; 
this rapid decay was understood to reflect the conversion of kinetic to potential energy 
under stable stratification, to which this scaling was assumed to apply exclusively. 2 

It has been suggested 3 that measurements 4 in very high Rayleigh number Rayleigh- 
Benard convection experiments are consistent with Bolgiano scaling. Although the evi- 
dence is not conclusive, this motivates the study of Bolgiano scaling in unstably stratified 
flow. L’vov and Falkovich have observed 5 that whereas Kolmogorov scaling corresponds to 
constant energy flux e, Bolgiano scaling corresponds to constant N, which can be identified 3 
with entropy flux. Their analysis also supports the suggestion 1 ’ 2 that Boussinesq turbu- 
lence typically exhibits both regimes, with Bolgiano scaling dominant at large scales and 
Kolmogorov scaling dominant at small scales. 

Bolgiano scaling will be considered here from the viewpoint of the Yakhot-Orszag 
(YO) theory. 6 It will be shown to arise from forcing of the temperature equation alone. 
A similar calculation for MHD has been made by Fournier et ai. 7 However, in this case 
forcing of the magnetic field equation leads to a Kolmogorov spectrum for the case of 
constant magnetic energy flux. This apparently occurs because the coupling constant is 
dimensionless in MHD but dimensional in Boussinesq turbulence. 

II. Langevin models for Boussinesq turbulence 

The Boussinesq equations in Fourier space are 

(-ZU? + vk 2 )ui = J dp u m (k - P)u n (p) + gP i3 (k)T(k) 

(-W + Kk 2 )T = -iiJfep J dp[u p (k - p)T(p) + u p (p)T(k - p)] (1) 
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In Eq. (1), the standard notation, 


k = (u>,k) = {u,ki,k 2 ,k 3 ) 

Pij( k) = Sij - kikj/k 2 

■P»mn(k) = &77i-P»n.(k) d" k n Pi m (k) 

is used. The buoyancy term coupling u and T is direction dependent. However, if Bolgiano 
scaling can arise independently of stratification, this dependence should be considered 
characteristic of the production mechanism only and as absent or irrelevant in the inertial 
range. Then it is not unreasonable to investigate an isotropic theory by replacing T in Eq. 
(1) by a convected vector field b{. Then 

( iw d - vk Jit, ~ 2 ®^*n»n(k) ^ dp u m (k p)u,i(p) d- P-fip(k)&p 

(-iw + K,k 2 )bi = ~^ik p J dp[u p (k - p)bi(p) + u p (p)bi(k - p)] (2) 

L’vov 3 suggests that the locality of inertial range interactions, expressed analytically as 
the finiteness of the DIA integral expression for N developed later in Sect. Ill, supports the 
possibility of Bolgiano scaling in unstably stratified flow: locality implies that inertial range 
scales interact among themselves and not with production range scales. This property of 
locality might also be invoked to justify Eq. (2) as a model of Boussinesq turbulence. 
Nevertheless, some caution in replacing the Boussinesq equations Eq. (1) by the isotropic 
model Eq. (2) is required, because vanishing of the correlation < u\T ' > in isotropic 
turbulence is a kinematic requirement for Eq. (1), whereas the correlation < u\b\ > in 
Eq. (2) need not vanish in this case. Accordingly, we leave more complete consideration of 
the relation between true Boussinesq turbulence and the isotropic model Eq. (2) for later 
investigations. 

The direct interaction approximation (DIA) for Boussinesq turbulence has been formu- 
lated by Kraichnan 8 . DIA for Eq. (2) can be written as a system of generalized Langevin 
equations 9 

ds 77(k,t,s)ui(k,s) d- /“( k,i) d-0&j(k,t) 

&i(k,t) = f ds ri b (k,t,s)bi(k,s) + f ds rj bu (k,t,s)ui(k,s) d- /f(k,<) (3) 

Jo Jo 
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where the damping functions axe defined by 

7?(k, t,s)=^ [ dp dq P rmn (k)P mr , (p)Pn,(q)G(p, t, s)Q u (9, t, s) 

4 dk=p+q 

77 6 (k,f,s) = j f dp dq k p qiP lp (p)G b (p,t,s)Q u (q,t,s ) 

^ «/k=p-fq 

77 5ri (k,t,s) = \ f dp dq { k p P pr9 (p)G(p,t,s)Q'l; r (q,t,s ) 

4 dk=p+q 

+ k p q r G b (q , i, (j>> «)} 


( 4 ) 


and the force correlations are 


< 


/“/“> = 1 / dp dq P imn (k)P ir4 (k)P mr (p)P„,(q)Q u (p,f,s)Q u (9,<,5) 

^ */k=p+q 

< /»/> > = 1 / dp dq Pw(k)fc p {P mp (p)Q tt (p,<,s)Q^(9,t,5) 

«/k=p+q 

+ P np (q)Q^(p,M)Q“(9,M)} 

< f b f b > = \ f dp dq {k p kiP lp (p)Q u {p,t,s)Q b (q,t,s) 

^ t/k=p+q 

+ 5 )} 

The scalar Green’s functions G and G b in Eq. (4) are defined as is usual in DIA by 

1 Sui 

fl 1 U 

a= d^Uff Pi ’ 

Qb _ 1 Mi - 

- difr 

In these equations, the correlation functions axe defined by 

< tii(<)uj(s) > = J dp Qij(p,t,s ) 

< Ui(t)bj(s ) > = J dp Qijip^s) 

< bi(t)bj(s) >=[ dp Q^(p,i,s) 


(5) 


J 


Isotropy implies that 


Qij(p>M) = Q u (p> t > s ) p ij( p) 
Qij( P,M) = Q b (p,i,s) 6 ij 
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The Y 0 theory can be understood as a solution of these equations in a special limiting 
case. First, the geometric coefficients axe all evaluated in the distant interaction limit 
k/p and k/q — > 0. Second, the triangle integrals axe all regularized by restricting the 
integration to the region in which either p > k or q > k. Finally, the time dependence 
is also evaluated in the distant interaction limit. Namely, assuming arbitrary scaling in 
the inertial range G(p, t,s) = G(ap r (t — s)) for some positive exponent r and suitable 
dimensional parameter a, the similarity form for G satisfies 

0 ( 0 ) = 1 

G(r) < 1 for t > 0 
G(r)dT < oo 

therefore 

p r G(p r (t - s)) ~ S(t — s) for p —* oo 

by standard properties 10 of delta functions. In Eqs. (4) and (5), this same argument 
is applied to products of Green’s functions and time correlation functions in the distant 
interaction limit. Therefore, in this limit, the damping is Markovian and the forcing is 
white noise in time. Fluctuation dissipation relations relating the correlation and response 
functions are valid in this limit. 

In Bolgiano scaling, Eq. (5) implies that the forcings satisfy the scaling relations 

< /“/“ > ~ < uA > ~ AT 19 / 5 

( 6 ) 

<f b f b > ~ k~ 3 

The velocity force spectrum is steeper than the k~ 3 force corresponding to a Kolmogorov 
spectrum; 5 accordingly, Kolmogorov scaling will apply at sufficiently small scales. The 
Bolgiano scale ks can then be understood as the scale beyond which Kolmogorov forcing 
dominates the velocity equation. Thus, 

ek i 3 - 5 6 / 5 iV 3 /5j bs 19/5 

Since in view of Eq. (6) the DIA velocity force scales like the cross correlation in 
Bolgiano scaling, it is reasonable to assume that the forcing of the velocity equation is 
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due entirely to the coupling to the b field. With these assumptions, DIA reduces to the 
” correspondence principle” of the Yakhot-Orszag theory 6 for Bolgiano scaling in Boussinesq 
turbulence, namely that the Bolgiano scaling regime is described by the noninear Langevin 
model 

( — H - uk )uj = -iPimn(k) U m {k - p)«n(p) + gPip(k)bp 

(— uo + nk 2 )bi = — J dp[u p (k — p)bi(p) + u p (p)bi(k — p)] — 7 k 2 Ui + fi ( 7 ) 

in which the force correlation is 

< > = D 8 (k + k'^ijk-v (8) 

and an arbitrary scaling exponent replaces the value y = d appropriate for Bolgiano 
scaling. 6 

The additional transport coefficient 7 corresponds to the damping function T] bu in 
the DIA Langevin model. This added term does not affect the RG calculations of v and 
k. The physical case, defined so that the stirring coefficient D has the same units as 
N, corresponds to y — d and e = 8. Although the damping terms have been replaced 
by the original nonlinear terms of the equations of motion, their effect will be evaluated 
perturbatively in the same ’’one loop” approximation used in the DIA Eq. ( 4 ). 

Suppose that this procedure were applied, not to the model equation ( 2 ) for the vector 
field b, but to the original equation (1) for T. Then it is easy to check that the corrections 
to viscosity and diffusivity would be tensors. However, the consistency of this calculation 
is uncertain. For as inspection of the DIA Langevin equations shows, anisotropic damping 
and correlations must correspond to an anisotropic force correlation with nonzero velocity 
and cross correlations. Only in the case of an isotropic theory do symmetry conditions 
dictate the form of the force correlation. But it is also easily verified that the isotropic 
part of the viscosity and diffusivity corrections computed in this manner would coincide 
with the results of the present calculation. 

The RG calculations follow as in Ref. 6. They lead to recurrences 

~ [A(a)(# - d+ 6 - e) + B(a) - C(a)] ( 9 ) 

S = ^ (1 ' (10) 
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in which 


e = 8 + y - d 


a — k/ v 


D = 2irDSd 

As in the Y O analysis of the passive scalar, the coefficients A, • • • , D axe rational functions 
of a: 


A(a) = 
B(a) = 

CM = 

D(a) = 


c - (-2 

4a(a + 1) 2 
c j i 2 -f- 3a -|” 4 
4a(a + l) 3 
7a 2 + 21a + 12 
4a(a + l) 3 
2a + 1 
4a 2 (a + l) 2 


( 11 ) 


Solving the recurrence relation for v when d = 3, 

- 4) = 5 D[A(a)(<P - i + 6 - «) + B{a) - C(«)p- - k~ B ') 

The initial conditions are v = vb when k = ks correspond to the highest wavenumber 
governed by Bolgiano scaling. The effective viscosity ug will be given by Kolmogorov 
scaling and not in general by the molecular viscosity. 


In the Y 0 theory, an expansion is made about the logarithmic theory e = 0 recognizing 
that evaluation of integrals beyond the one loop level will introduce further corrections of 
order e n with to > 1 to the constants. At lowest order, consistency demands that all 
such corrections be dropped, including those already computed in the present one loop 
approximation. Dropping these corrections and setting 6 = 8, 

where 

A = A(a)(d 2 -d + 6) + B(a) - C(a ) 


The quantity a, an inverse Prandtl number, is evaluated by assuming it is scale independent 
and can be evaluated by setting the viscosity v proportional to the conductivity k in the 
inertial range. This leads to a quartic equation for a 

3 a 4 + 9a 3 - 2a 2 - 15a - 5 = 0 (12) 
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which has only one positive root 


a = 1.291 


(13) 


Then from Eq., 

A = .7929 

It must be stressed that in this procedure, e is never set to any value other than 
8. These results could also be derived directly from the simplified DIA Langevin model 
without introducing recurrence relations following Woodruff. 11 By regularizing the triangle 
integrals, crossover to a do minant sweeping interaction 11 when e = 7 is suppressed. 

The inverse Prandtl number of Eq. (13) is very dose to the value 1.3929 obtained 
for forced convection by Yakhot and Orszag 6 (compare also Ref. 7). The close proximity 
of these n umb ers is consistent with a common assumption 12 that the turbulent Prandtl 
number is independent of Richardson number in the atmospheric surface layer. Although 
experimental evidence does not support this assumption unambiguously, 12 the theoretical 
value of Eq. (13) is within the range 1.00 < a < 1.43 observedin some recent measurements 
of buoyant plumes. 13 

Spectra are defined as usual by 


E.(p) = \ j 

f dp 

r & oh® 

J — oo 

Eh(p) = \ j 

f dp 

/ du QUp) 

J — oo 

E i (p)=l j 

f> dp 

r d»Q i u (p) 

J — oo 


and evaluated at lowest order in perturbation theory as 

E, = ^ I dp r du I G(p) I 2 1 G t (p) I 2 g 2 Dp~* 

= id p[°°du G(p) | Gi(p) I 2 gDp~< 

- ^(5A)- 2 /V/‘d 2/ V’ /2 
= i(?A)-GV 2 / ! D‘/V ,/s 
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Then 


u(k) = a^D^k^ 5 

E u (k) = a u g 4 / 5 D 2 / 5 *" 11 / 5 

E h (k) = o^/ 5 D 3 / 5 r 9 / 5 (14) 

E b (k) = a^- 2 / 5 D 4 / 5 fc- 7 / 5 

where 

a v = 0.440 
a u = 1.310 
a h = 0.576 
a b — 0.679 

Introducing an inverse integral scale kf so that the energy is contained entirely in the 
inertial range kf < k < oo, the single point turbulence quantities axe 

K = i < u p u r > = 5„„ s 4 /5 D v= t ;'/ 5 

E = ^ < u p b p > = 5 ot9 i/5 D V5j.J 4/5 (15) 

c = i > = ^ w - 2/5 d‘/ 5 i; 2/s 

= v{kf) — o. v g x ^ H 2 / 5 kJ 6 / 5 

In these equations, corrections due to a Kolmogorov regime at small scales have been 
ignored. 

III. Elimination of the force amplitude 

The expressions for spectra and single point turbulence quantities developed above 
all contain the quantity D. In the YO theory, this quantity is expressed in terms of the 
dissipation rate N through the DIA inertial range energy balance, 14 

0 = ik p / dpdq < [u p (p,t)6i(q,t) +u p (q,t)5i(p,t)]6j(-k,t) > 

*/k=p+q 
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Assuming a steady state, let 


S^k^q) = ik p < [u p (p,t)fe»(q,<) + « p (q,*)6i(p,<)]&i(-k,i) > 


Integrating the DIA energy balance over all time separations, 

S(k,p, q) = 2Q{k,p,q)k p P pmn {p)P m i{^)Pni{qdQ h { k )Q h ( < l ) 

+ 0(fc,p,9)[Jfc P ? f ^ i (k)Ppi(p)Q' l (^)Q' l (p) - k p k l P pi (p)Pu(q)Q h (p)Q h (<l)] (16) 

+ Q{k,p,q)[k p qiPi P (p)6u{k)Q u {p)Q b {k ) - k p k l P pl (p)S ii (q L )Q u {p)Q b {q)} } 

+ (?> ?) 


where (p, q) denotes interchange of wavenumbers p and q and the time scale 0 is evaluated 
using fluctuation dissipation relations, valid in the present case of Markovian damping and 
white noise forcing, 

<Ui(t)bj(s)> = J dsG(k,t,s) < Ui(s)bj(s) > 

< Ui(s)bj(i) > = J dsG b (k,t,s) <Ui(s)bj(s) > 
for t > s. The result is 

®{k,p,q)= f da G u (p,a)G b (k,a)G b (q,a) 

Jo 


Following Kraichnan, 14 

N=[ dk f dp dq S(k,p,q) - / dk f dp dq S(k,p,q) (17) 

Jk<k 0 Jp,q>ko J2k 0 >k>k 0 Jp,q<k o 

must be finite and independent of ko when inertial range scalings are substituted in 5 . The 
triangle condition k = p + q is understood in the integrals. Convergence must be checked 
in the cases p —* 0 and q — ► 0. We find 

lim,_ 0 £ d<\ S(k,p,q) = 0(q 4 )Q{q) 


therefore 


/ S(k,p,q)~fdq 
J q* s/0 v 


2-ol 


< oo 
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because in all integrals, a < 3. 

To evaluate the flux integral, define the quantity T(k,p,q ) by 

/ dp dq T(k,p,q) = / dp dq S(k,p,q) 
J */k=p+q 

The sc aling s are such that 


T(\k,\p,\q) = \- 3 T(k,p,q) 

a result equivalent to the choice of a k~ 3 force. This implies 15 that the inertial range trans- 
fer does not depend on the wavenumber ko and permits the transfer integral to be reduced 
to a double integral over a finite region following Kraichnan. 14 Numerical integration of 
Eqs. (16)-(17) leads to 

N = 4.01D 

This result can be used following YO to eliminate the force amplitude D from Eq. 
(14) for the spectra and Eq. (15) for single point quantities. The results for the scale 
dependent viscosity and spectra are 

u(k) = O.Z3Zg 2/5 N 1/5 k- i/5 
E u (k) = 0.752 g 4 / 5 N 2 / 5 k~ n / 5 

E h (k) = 0.250$ 1/5 1V 3 / S fc~ 9/5 (18) 

E b (k) = 0.336 g- 2 / 5 N 4/5 k~ 7 / 5 

For the single point turbulence quantities: 

K = 0.627g 4 ' 5 N 2 / 5 kJ* /5 

H = O.3130 1 / 5 JV 3 / 5 fc J 4/5 (19) 

C = 0M0g~ 2/5 N 4 / 5 kJ 2/s 

These results lead to a formula for the turbulent viscosity in terms of single point 
quantities, 

KO 

v T = 0.633— (20) 
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which can be applied in transport models. In problems in which the heat transfer H is 
constant, the relations 

H 2 = 0.18 6KC 

N = 12.8 gK~ 2 H 3 = 1.03 gC 3/2 K~^ 2 

may also be useful. 

IV. The transport coefficient 7 

The transport coefficient 7 can be computed after v and k because it does not enter 
the recurrence relations for these quantities. We find 

d 1 1 gDp d ~ l 1 f 7a 2 + 4o + 1 
dp ~ 4 v 3 p 6 +y d(d+2y 4a 2 (a + l ) 3 

+ TTT~TtA3 f 7a2 + 4a + 1 - ( 3a + l)(a + 1)(6 + 2 /)] 

4a 2 (a + 1)'* 

, * 3a + l 1 2 

4a 2 (a + l) 2a(a + l) 2 a(a + l) 3 

Note that 7 does not scale like v and /c; instead, using v ~ p~ e / 5 , we find 

7 ~ p3c/5— (e— 2 ) 

indicating that 7 is irrelevant for e < 5. Accordingly, to evaluate the amplitude for 7 , it is 
appropriate to consider u and a as known, but to expand otherwise about e = 5, or y — 0. 
This procedure leads to 

7 = 0.00263ff -1/5 JV 2 / 5 Jfe~ 6/5 

It appears that this quantity can be neglected for practical purposes. 

V. Iterative solution of the DIA Langevin equations 

Let us compare the present calculation using the YO theory with the original DIA 
Langevin equation. The YO ’’correspondence principle” can readily be understood in terms 
of this Langevin model. The necessity of a stirring force is immediate from DIA, and the 
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k~ 3 scaling is equivalent to Orszag’s transparency condition 15 on the inertial range. The 
YO theory begins by generalizing this model: the DIA response equation is satisfied by 
a family of scaling solutions depending on the force scaling exponent y. Then the theory 
computes perturbatively about the logarithmic (Gaussian, asymptotically free) theory in 
which e = 8 + y — d= 0. DIA does not suggest this procedure since only e = 8 is consistent 
with the DIA response equation. However, if it is agreed that this procedure is appropriate, 
the calculation of amplitudes in the physical case is consistent: in this calculation, e is never 
set to any value other than 8. The numerical accuracy of this procedure remains uncertain. 


The stirring force amplitude is eliminated by appeal to the full DIA response equation 
which is neither regularized nor evaluated in the distant interaction limit. However, even 
this calculation is incomplete because the time dependence has been calculated assuming 
Markovian damping acting against a white noise in time forcing. Finally, the viewpoint 
of DIA, the correspondence principle for Boussinesq turbulence was incomplete because 
it assumed that f u could be neglected in comparison to < ub >. The consistency of this 
assumption has not been demonstrated. 


All of these considerations suggest that we have only computed a first approximation 
to DIA. In abbreviated notation, this solution is 


Gi(k,<r) ~ exp (— <rfc 2 / 5 ) 

<?$(*, <r) ~ exp(-<rib 2 / 5 ) 

Q“(jfe,<r) ~ k-^G^a) 

Qi(k,(r) ~ k~ 19 / 5 Gi(k,cr) for a > 0 
Qi{k,a) ~ k~ 19 ^G\{k,a) for a < 0 


where time stationarity has been assumed and <j = t — s. The DIA Langevin model can 
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be used to generate corrected damping and forcing: 



Vi{k,‘ 

cr) ~ 

/ GiQ * 


Vi(k, 

<t) ~ 

/<?!<?? 



<r) ~ 

/ G1Q1 

F uu 

=< fifi 

> ~ 

Q1Q1 

F ub 

V 

II 

> ~ 

Q 

pbb 

II 

A 

> ~ 

<??<?; 


With these corrected dampings and forcings, new response functions and correlation func- 
tions G2, <?2, can be calculated assuming that the damping and forcing are 

known: 

(?2 + 7]l * = 0 

G b 2 +v b i*G b 2 +v b i u *G b 2 =0 

QZ+V 1 *QZ = G 1 *Fr 

Q b 2 + vl*Q h 2 + rf u *QZ = G\*F? 

This iterative procedure answers the objections raised above. At the second itera- 
tion, the distant interaction limit will be corrected, the Markovian time dependence will be 
modified (hence the fluctuation dissipation relation will not be satisfied) and the complete 
DIA forcing will be computed. Regularization of the response equation cannot be avoided 
in this kind of theory. This procedure introduces an arbitrary parameter into the theory 
since we can also regularize by requiring p > ak for any positive constant a . 16 We have 
followed YO by setting a = 1 . 

APPENDIX: Integrals in calculation of renormalized transport coefficients 
The following notation is used: 

G(k) = (—iu> + uk 2 )~ l 
Gb(k ) = (—ia? + k fc 2 ) -1 
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Viscosity corrections axe calculated from 


J du £ dp 

(-^iG(k - p)P m „(k - p)g 2 P 3n (p)G{p)G(-p)G b (p)G b (-p)p~ y 

iG{p)P nr9 (p)g 2 P m8 (k-p)G(p-k)G(k-p)G b (p-k)G b (k-p ) | k-p | _v 

where f dp denotes integration over a d — 1 dimensional sphere of radius p. The integrals 
are evaluated in the limit k/p — * 0. The frequency integrations in this limit axe 

/ oo o— i 

- P)G(p)G(-p)G<,(p)G„(-p) = x 

•[(a + 2 )(o + 1) + k v p v p 2 (a 2 + 3a + 2)} 

and the second frequency integral can be found by replacing p by k — p in this result. 
Similarly, 

k 2 ^ = — -k p [ du> (f dp 
dp 4 J _ 00 J 

| G(p -k) | 2 | G b (p - k) | 2 G b (p)p m DS mp | k-p |-v 
I G(p) \ 2 \ G b (p) | 2 G b (k - p)Dp~ y S pm 

The frequency integrals axe evaluated by interchanging v and k, or equivalently by replacing 
a by 1/a. Also, 

k 2 ^- — —~kp f du> <£ dp 
dp 4 ' */_<*, J 

G(p - £) | G b (p - fc) I 2 G b (p)gDp m P P i ( k - p) | k - p \~ y 
+ G(p) | G b (p) | 2 G b (k - p)gD(k m - p m )Ppi(p)p~ v 
+ G(k - p)G(-p) | G b (p) | 2 gDP pmn ( k - p)P ni (p)p~* 

+ G(p)G(p - k) | G(p - k) | 2 gDP prnn (p)P n i(p — k) | p — k [“ y 

In evaluating products of projection operators in the distant interaction limit, the first 
order Taylor series 

Pij(p - k) = Pij(p) + p~ 2 k v P vij ( p) + • • • 

is useful. 


15 



REFERENCES 


X. R, Bolgiano, ’’Turbulence spectra in a stably stratified atmosphere,” J. Geophys. Res. 
64 2226 (1959). 

2. A. S. Monin and A. M. Yaglom, Statistical Fluid Mechanics vol. II, MIT Press, Cam- 
bridge (1975). 

3. V. S. L’vov, "Spectra of velocity and temperature with constant entropy flux of fully 
developed free convective turbulence,” Phys. Rev. Lett. 67 , 687 (1991). 

4. X.-Z. Wu, L. Kadanoff, A. Libchaber, and M. Sano, "Frequency power spectrum of 
temperature fluctuations in free convection,” Phys. Rev. Lett. 64, 2140 (1990). 

5. V. S. L’vov and G. E. Faikovich, ’’Conservation laws and two-flux spectra of hydro- 
dynamic convective turbulence,” Physica D 57, 82 (1992). 

6. V. Yakh ot and S. A. Orszag, "Renormalization group analysis of turbulence,” J. Sci. 
Comput. 1, 3 (1986). 

7. J. D. Fournier, P. L. Sulem and A. Pouquet, ’’Infrared properties of forced magneto- 
hydrodynamic turbulence,” J. Phys. A 15, 1393 (1982). 

8. R. H. Kraichnan, ’’The direct interaction approximation for shear and thermally driven 
turbulence,” Phys. Fluids 7 , 1048 (1964). 

9. R. H. Kraichnan, "Eddy viscosity in two and three dimensions,” J. Atmos. Sci. 33, 
1521 (1976). 

10. I. M. G elfan d and G. Shilov, Generalized Functions vol. 1, Academic Press, New York 
(1964). 

11. S. L. Woodr uff , ’’Dyson equation analysis of inertial range turbulence,” Phys. Fluids 
A 4, 1077 (1992). 

12. A. S. Monin and A. M. Yaglom, Statistical Fluid Mechanics vol. I, MIT Press, Cam- 
bridge (1975). 

13. A. Shabbir and W. K. George, ’’Experiments on a round turbulent buoyant plume,” 
NASA technical memorandum 105955 (1992). 

14. R. H. Kraichnan, ’’Inertial range transfer in two and three dimensional turbulence,” 
J. Fluid Mech., 47 , 525 (1971). 

15. S. A. Orszag, Statistical theory of turbulence in Fluid Dynamics, ed. R. Balian and 
J.-L. Preube, Gordon and Breach, New York (1977). 


16 


16. R. H. Kraichnan, ”An interpretation of the Yakhot-Orszag theory,” Phys. Fluids 30 
2400 (1987). 


17 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reoortinq burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing dal^urces, 
oathering^and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other asp^ of this 
collection of information including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information 0 perations and Reports, 12 5 Je 
Davis Highway, Suite 1204, Ariingion^VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1 . AGENCY USE ONLY ( Leave blank) * 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

May 1994 Technical Memorandum 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 


Renormalization Group Theory of Bolgiano Scaling in Boussinesq Turbulence 


6. AUTHOR(S) 


WU-505-90-5K 


Robert Rubinstein 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8877 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORINGMONITORING 
AGENCY REPORT NUMBER 

NASA TM-1 06602 

ICOMP-94-8 

CMOTT-94-2 


11. SUPPLEMENTARY NOTES , 

Robert Rubinstein, Institute for Computational Mechanics in Propulsion and Center for Modeling ot Turbulence and 

Transition, NASA Lewis Research Center (work funded by NASA Cooperative Agreement NCC3-233). ICOMP 
Program Director, Louis A. Povinelli, organization code 2600, (216) 433-5818. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b- DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words) 


Bolgiano scaling in Boussinesq turbulence is analyzed using the Yakhot-Orszag renormalization group. For this 
purpose, an isotropic model is introduced. Scaling exponents are calculated by forcing the temperature equation so that 
the temperature variance flux is constant in the inertial range. Universal amplitudes associated with the scaling laws are 
computed by expanding about a logarithmic theory. Connections between this formalism and the direct interaction 
approximation are discussed. It is suggested that the Yakhot-Orszag theory yields a lowest order approximate solution 
of a regularized direct interaction approximation which can be corrected by a simple iterative procedure. 


14. SUBJECT TERMS 

Buoyancy; Free convection; Turbulence; Renormalization group 

17 SECURITY CLASSIFICATION I 1 8. SECURITY CLASSIFICATION 1 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 

NSN 7540-01-280-5500 


15. NUMBER OF PAGES 

19 


16. PRICE CODE 

A03 


, 20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 







